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Common beans (Phaseolus vulgaris L.), like many legumes, are rich in iron, zinc, and 
certain other microelements that are generally found to be in low concentrations in 
cereals, other seed crops, and root or tubers and therefore are good candidates for 
biofortification. But a quandary exists in common bean biofortification: namely that the 
distribution of iron has been found to be variable between the principal parts of seed; 
namely the cotyledonary tissue, embryo axis and seed coat. The seed coat represents ten 
or more percent of the seed weight and must be considered specifically as it accumulates 
much of the anti-nutrients such as tannins that effect mineral bioavailability. Meanwhile 
the cotyledons accumulate starch and phosphorus in the form of phytates. The goal of 
this study was to evaluate a population of progeny derived from an advanced backcross of 
a wild bean and a cultivated Andean bean for seed coat versus cotyledonary minerals to 
identify variability and predict inheritance of the minerals. We used wild common beans 
because of their higher seed mineral concentration compared to cultivars and greater 
proportion of seed coat to total seed weight. Results showed the most important gene 
for seed coat iron was on linkage group B04 but also identified other QTL for seed coat 
and cotyledonary iron and zinc on other linkage groups, including B11 which has been 
important in studies of whole seed. The importance of these results in terms of physiology, 
candidate genes and plant breeding are discussed. 



Keywords: advanced backcross breeding method, cotyledon, embryo axis, iron and zinc concentration, use of wild 
beans, seed coat 



INTRODUCTION 

Biofortification is a relatively recent addition to breeding goals in 
plants based on improving the nutritional quality of the edible 
portion of the plant through traditional or transgenic approaches 
(Dwivedi et al., 2012). To date, most biofortification work has 
concentrated on micronutrients and vitamins, although con- 
ceivably protein content, amino acid distribution and beneficial 
secondary metabolites could be considered to be goals of biofor- 
tification (Welch, 1999). The stated goals of most biofortification 
work is to reduce mineral or vitamin deficiencies and where 
needed protein deficiencies. These deficiencies manifest them- 
selves in conditions of iron deficiency anemia (IDA) but are not in 
themselves actual diseases but rather imbalances in the diet that 
need to be approached through modifications in the diet (Pfeiffer 
and McClafferty, 2007). 

Micronutrient deficiencies are predicted to affect half the 
world's human population, with IDA being an especially com- 
mon health concern affecting at least 2 billion people. IDA is 
caused by low consumption of iron especially in reproductive 
age women and developing adolescents (Welch, 1999). Zinc defi- 
ciency is suspected to be equally as common but has not been as 
well documented as IDA (Welch and Graham, 2002). While IDA 



causes losses in work productivity and developmental problems, 
zinc deficiency causes lowered disease immunity and stunting. 
These types of deficiencies are sometimes difficult to address 
through supplementation or fortification for technical, societal 
or economic reasons and therefore these minerals do not always 
reach the poor consumer (Dwivedi et al., 2012). How to reach 
the bottom of the pyramid of societies' economic strata (people 
making less than 2 dollars a day) with food containing sufficient 
micronutrients has been debated but one way is through biofor- 
tification of staple crops that are consumed in high amounts by 
the poor, such as rice, wheat and beans (seed crops) or potato and 
cassava (root and tuber crops). 

Common bean is a highly nutritious food because of its 
balance of carbohydrates to protein (between 4:1 and 3:1), con- 
tent of important vitamins and its micronutrient concentra- 
tion, especially for iron which is at much higher levels than 
in the starchy staples of barley, corn or wheat. Common bean 
is widely grown in many parts of Asia, Africa, Europe and 
North, Central and South America and has two major types, 
the large seed Andean beans and the small seed Mesoamerican 
beans (Broughton et al., 2003). Consumption of common beans 
generally is higher than 20 to 30 kg/year in rural northeast 
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Brazil, Central America and Mexico and can be as high as 40 
to 60 kg/year in regions where meat is scarce such as in the 
Great Lakes of Africa. In Europe and North America consump- 
tion is low at 5 kg/year or less. Like other legumes, common 
bean is usually combined in the diet with a starch based food 
either in a mixture of two components (pulse + grain) or as 
a side dish such as "dais" made of lentils or "Frijoladas" made 
of beans. In some places common bean is the national dish as 
in Brazil where "Feijoada" is served daily or in the Dominican 
Republic where the "Bandera" is a twice-a-day meal made of 
rice with a side of beans. Certain regions of Central and South 
America consume mashed or re-heated "Calentao" beans as 
breakfast. 

One reason for the micronutrient density of legumes is their 
anatomical seed structure, where a thick, maternally-derived 
seed coat or testa surrounds the expanded cotyledons of a 
quiescent but fairly large and well developed embryo. This 
structure differs radically from cereals where a thin, maternally- 
derived aleurone layer surrounds specialized endosperm tis- 
sue with a less well developed embryo in the bran. Breeding 
programs for cereals have had little concern for the sub- 
compartmentalization of micronutrients given that the bran is 
removed in milling processes and the major target tissue is 
the endosperm. However, in legumes, for the most part, the 
seed is consumed whole after a process of boiling, rather than 
by any processes of milling or grinding. In addition, the seed 
coat makes up 7 to 10% or more of the total seed weight of 
beans and is a source of consumer preference, the cotyledon 
makes up 85% or more of total seed weight and the embryo 
is only 2 or 3% of seed weight but is dense in nutrients 
(Ariza-Nieto et al., 2007). 

Seed size varies more than seed coat thickness so smaller 
seeded legumes especially wild relatives of beans tend to have even 
higher percentages of seed coat compared to cotyledonary tissue 
of cultivated beans. The presence of these different seed tissues 
means that all three components of the common bean seed should 
be targets of biofortification: seed coat, cotyledons and embryo. 
An understanding of mineral distribution, loading and inheri- 
tance into each tissue is essential for making progress in breeding 
of this crop. 

The goal of this research, therefore, was to evaluate the concen- 
tration of iron and zinc, in particular, and other minerals more 
generally in seed coats that were separated from cotyledons for 
a common bean population derived from a wild (small-seeded) 
common bean crossed with a cultivated (large-seeded) common 
bean. This analysis was part of a biofortification breeding pro- 
gram using the advanced backcross (AB) quantitative trait loci 
(QTL) approach as outlined in Tanksley and Nelson (1996); Blair 
et al. (2006) and Blair and Izquierdo (2012). The identification 
of QTL for seed coat versus cotyledonary mineral concentrations 
can be used in molecular breeding of common bean as regular 
peeling of seed coats for analysis is time consuming and onerous 
for the researcher. In addition, we were interested in determining 
the value of wild beans as a source of higher iron concentration 
as they have been found to have very high levels of this min- 
eral (Guzman-Maldonado et al., 2004) but have a disadvantage 
of small seed size. 



MATERIALS AND METHODS 
SEED SOURCE AND FIELD CONDITIONS 

An advanced backcross (AB) population of 138 BC2F2:5 lines 
(cross NH 21154_ (F2)— 1— M— M) was developed by crossing the 
cultivated recurrent parent Cerinza with the wild donor parent 
G10022, a source of high seed mineral concentration as described 
in Blair and Izquierdo (2012). Cerinza as the recurrent or recipi- 
ent parent is large-seeded and of red color, and is rounder than a 
kidney bean. 

Cerinza, like other Andean large-seeded beans, weighing 
approximately 50 g per 100 seed, has a good baseline of iron con- 
centration averaging 60ppm in multiple field trials (Astudillo 
and Blair, 2008). Meanwhile, the wild donor parent G10022 from 
Mexico is small seeded, weighing approximately 5 gper 100 seed 
and was discovered in a screening program for nutritional traits 
in wild beans and was found to have a very high concentration 
of 110 ppm iron concentration (Islam et al, 2002). Zinc concen- 
tration was less contrasting in previous studies of the two parents 
with G10022 having 38 ppm and Cerinza having 25 to 27 ppm. 

The advantage of using the AB population instead of a simple 
cross between wild and cultivated beans for the nutritional anal- 
ysis was that most of the genotypes were of similar seed size and 
adaptation (Blair et al., 2006; Blair and Izquierdo, 2012). They 
could therefore be grown together in the same field experiment at 
a site in Darien, Colombia (3°54'N, 76°30'W, 1485 m above sea 
level, average yearly temperature 20° C, average relative humid- 
ity 80%, average yearly rainfall 1288 mm) with moderately acid, 
loam soils (pH 5.6, Andisol). The field experiment (a random- 
ized complete block design with three replicates and 3 m long 
single-row plots) was managed with two foliar applications of zinc 
and boron as microelements (300 g ha^^ as chelates) carried out 
at 14 and 21 days after planting. Harvesting was by hand into 
clean plastic buckets to avoid soil contamination. Foliar disease 
pathogens were controlled with the fungicides, mancozeb, and 
benomyl, at doses of 1.0 and 0.3 kg ha^^, respectively. 

A sample of the harvest from the first two replicates was care- 
fully placed separately into two paper bags to keep them clean 
through transport. Back in seed laboratory, the seeds were hand 
washed with 70% ethanol to remove dust and dried in a station- 
ary oven. The seed were all weighed upon reaching 10% moisture 
and the seed color was recorded, although most of the genotypes 
had medium to large red seed (45 to 55 g per 100 seed) like the 
recurrent parent Cerinza. The similarity of the lines was a result 
of backcrosses of the recurrent cultivated parent as a recipient of 
the wild donor parent genes. 

MINERAL EVALUATION AND DATA ANALYSIS 

For each repetition, 12 g of seed were washed with sterile double- 
distilled water and peeled by hand using a sterile scalpel to 
remove the seed coat from the cotyledons. The separate seed 
coat and cotyledon samples were placed in separate envelopes 
and then dried for 24 h at 45°C in a bench-top oven, before 
grinding in a modified Retsh mill with 24 sample slots using 
zirconium grinding balls and Teflon grinding chambers. Seed 
coat and cotyledonary tissues were dried and weighed out into 
two replicates of ~0.25 g dry weight each which were analyzed 
separately before averaging. These samples were analyzed with 
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nitric acid-perchloric acid digestion as described in Blair et al. 
(2009). 

Digested samples were taken to dryness and re-suspended in 
15 mL of trace-metal grade nitric acid (2% v/v), prior to multi- 
element analysis by inductively coupled plasma — optical emis- 
sion spectroscopy (ICP-OES) (CIROS ICP Model FCE12; Spectro, 
Kleve, Germany) at the USDA-ARS Children's Nutrition Research 
Center at the Baylor College of Medicine (Houston TX). The 
instrument was calibrated daily with certified standards. Tomato 
leaf standards (SRM 1573A; National Institute of Standards and 
Technology, Gaithersburg, MD) were digested and analyzed along 
with the seed samples to ensure accuracy of the instrument's cali- 
bration. Because of the large amount of work involved in peeling 
sufficient seed, only 60 of the full set of 138 AB representing a 
random selection of lines for full-seed iron concentration were 
evaluated. These included 20 high, 20 moderate and 20 low seed 
iron containing lines for the cross covering the full range from 
60 to 95ppm (Table SI) as measured for Darien and Palmira by 
Blair and Izquierdo (2012). In addition, the two parents, Cerinza 
and G 10022 were used as control genotypes twice in the analy- 
sis. Because of the higher seed coat ratio in G10022 only 4.5 g 
of this wild bean was needed to obtain sufficient seed coat for 
analysis. 

All quantitative data were analyzed using a general linear 
model and an analysis of variance for a split plot design for 
genotypes within tissue types analyzed in the software pack- 
age Statistix v. 8.0 (Analytical Software Inc.) and means were 
estimated to use for subsequent analysis. Population distribu- 
tions were evaluated for normality using the same software. 
Transformations were carried out with natural logarithm in cases 
were population results were skewed before use for QTL anal- 
ysis. Pearson's correlations were estimated based on the non- 
transformed data. QTL analysis was carried out with the means 
described above and the genetic map built by Blair and Izquierdo 
(2012). QTL were identified using two software programs QTL 
Cartographer v. 2.5 (Basten et al., 2001) for composite interval 
mapping analysis (CIM) and MapDisto v. 1.7 (http://mapdisto. 
free.fr/) for single point analysis (SPA). In the CIM analysis we 
performed the analysis every 1 cM (walkspeed) with a window 
size of 10 cM and using ten background markers in a forward- 
backward stepwise multiple linear regression model. In terms 
of population type, the B12 genetic model was assumed for 
the CIM analysis. Meanwhile a simple regression model was 
assumed for the SPA analysis. In both cases the homozygous 
donor parent allele class was combined with the heterozygous 
genotypic class. Significance thresholds were set at LOD 3.0 
{P < 0.001). 

The phenotypic variance controlled by a given QTL was deter- 
mined by its determination coefficient (R2), as defined by the 
software program. QTL for micronutrient concentrations were 
named using the mineral name, the abbreviation for seed coat or 
cotyledon and a two number code derived from the linkage group 
and the number of the QTL identified on that linkage group, sep- 
arated by a period. Genetic maps and QTL locations were drawn 
with MapChart v. 1.0 software program (Voorips, 2002) where 
map distances were reported in centiMorgans (cM) estimated 
with the Kosambi mapping function. 



RESULTS 

MINERAL VARIABILITY IN SEED COAT vs. COTYLEDONS 

A total of ten minerals were analyzed by the ICP-OES experi- 
ment for the two subsamples namely the seed coat and cotyledon 
(Table 1). These minerals were boron (B), calcium (Ca), cop- 
per (Cu), iron (Fe), potassium (K), magnesium (Mg), manganese 
(Mn), phosphorus (P), sulfur (S), zinc (Zn). The concentrations 
were calculated in part per million (ppm) equivalent to g/kg or 
(jLg/g. The seed coat was found to have low concentrations com- 
pared to the cotyledon for B, Cu, K, Mn and especially P and S. 
Meanwhile the seed coat had higher concentrations of Ca, Fe, Mg, 
and Zn. 

For the minerals of greatest interest to our study, average Fe 
and Zn concentrations of the cotyledon were similar to that of 
seed coat across the entire population (Table 1). However, the 
range of Fe in seed coat was very large with a minimum of 20 ppm 
to a maximum of 263 ppm, while Zn in seed coat ranged from 
17 to maximum of 54.5 ppm. By comparison cotyledonary val- 
ues were less variable ranging from 54.5 to 93 ppm for Fe and 
27 to 41 ppm for Zn. For other nutritionally important minerals, 
Ca ranged from 12500 to 21500 ppm in seed coat which was sur- 
prising given average of 284 ppm in cotyledon. The opposite was 
true for P which averaged 4332 in cotyledon but only 424 ppm in 
seed coat. Similarly S averaged 1982 ppm in cotyledons but only 
350 ppm in seed coat. The differences in the concentration aver- 
ages between seed coat and cotyledon for each of these minerals 
were significant in each case. 

Observation of population distribution types for iron in 
cotyledons and zinc in both tissues showed normality and a lack 
of high kurtosis or skewing (Figure 1, Table 1). The exception to 
this was the binomial distribution and high K and S values for 
seed coat iron suggesting simpler inheritance and perhaps a major 
QTL for this trait. This was in contrast to the normal population 
distributions for zinc concentrations in seed coat or for iron con- 
centration in the cotyledons. Therefore the inheritance of these 
minerals in these tissues appears to be controlled by multigenic 
or quantitative. The other minerals were also graphed for their 
population distribution (data not shown). The range, median, 
means and K or S values all indicate normal distributions and 
quantitative inheritance for the most part, although the -K" and S 
values were high for seed coat Cu and for cotyledonary B and Mg. 
Kurtosis was also observed for cotyledonary S (Table 1). 

STATISTICAL ASSOCIATIONS BETWEEN MINERALS 

In the ANOVA results, genotype effects for the advanced back- 
cross lines were highly significant {P = 0.01-0.001) in each case 
except B, for which the effect was lower {P = 0.05) while tis- 
sue effects were either highly or moderately significant (Table 2). 
Tissue X genotype effects were also highly significant for all min- 
erals. Replication differences were minor as evidenced by low 
sums of squares values (data not shown), which is explained by 
the good repeatability of the study, with separate samples man- 
aged from the field to the seed room to lab analysis with very 
similar treatment. Coefficients of variation (CVs) were below 3%, 
except for B and Ca which had higher CVs for the two error terms 
used Error (Replication x Tissue) to test significance of Tissue 
differences and Error (Replication x Tissue x Genotype) to 
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Table 1 | Descriptive statistics for seed coat vs. cotyledonary minerals evaluated as measured by inductively coupled plasma (ICP) analysis for 
the advanced backcross population of common beans. 
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' Mineral abbreviations: B, boron; Ca, calcium; Cu, copper; Fe, iron; K, potassium; IVIg, magnesium; Mn, manganese; phosphorus; S, sulfur; Zn, zinc. 

Tissue Abbreviations: coat, seed coat; coty, cotyledon (and embryo). Tlie minerals of most interest, iron and zinc are indicated in the columns which are in bold text. 

^Mineral concentrations given in part per million (ppm); note difference between micronutrients (B, Cu, Fe, Mn, Zn) vs. other minerals (Ca, K, Mg, P and S). 



test significance of Genotype and Tissue x Genotype difl^erences 
(Table 2). 

Parental comparisons showed that for the most part tissue dif- 
ferences were highly to moderately significant as were genotype 
differences, except in the case of the first mineral, B (P > 0.05), 
which was non-significant for any of the sources of variation 
(Table 2). Tissue x genotype effects were also non- significant for 
Cu and Zn amongst the parents but were significant for all other 
minerals although to a lesser extent and lower probability than for 
the advanced backcross lines. 

Similarly, probability values for differences between the tis- 
sues were less significant than for the lines in most cases although 
this was more notable for the minerals B and Fe among the lines 
(differences at P < 0.05) with all other minerals surpassing the 
level of high significance (differences at P < 0.001). Some of the 
genotype comparisons for minerals were very highly significant 
(differences at P < 0.0000). Tissue x genotype effects were less 
observable for B but were also of very high significance for the 
remainder of the minerals (differences atP < 0.0000). 

The average parental values for Fe and Zn (Figure 1, arrows) 
were of interest given our attention to these micronutrients. 
Most noticeably, the wild donor parent (G10022) had more than 
double the seed coat Fe concentration (93.28 ppm) of the cul- 
tivated recurrent parent (Cerinza, 40.52 ppm). Meanwhile the 
cotyledonary Fe concentrations were similar but reversed with 



Cerinza having higher iron in the cotyledons (80.16 ppm) than 
the cotyledons of G10022 (67.51 ppm). Differences were signif- 
icant based on paired t-tests, but higher for seed coats (P = 
0.0039) than for cotyledons (P = 0.0180). For Zn concentrations 
the cotyledons were about 6-8 ppm richer in this mineral than 
the seed coats. However, in both cases G 10022 had more Zn 
(28.7 ppm in seed coat, 36.2 ppm in cotyledons) than Cerinza 
(26.5ppm in seed coat, 32.16ppm in cotyledons). Paired f-tests 
showed high significance for the contrast in seed coat (P = 
0.0033) but non-significance for the contrast in the cotyledon 
(P = 0.0684). 

The Pearson's correlation coefficient results (Table 3) con- 
firmed tissue X genotype interaction and that the relationships 
among minerals was different in the two types of tissues, the seed 
coat and the cotyledon. For cotyledonary tissue, Mn was corre- 
lated with P and S just as K was correlated with Mn but none 
of the negative correlations were highly significant. Less signifi- 
cant (P < 0.05) and negative correlations of Ca and K or Fe and 
Mg were observed in the cotyledons. Likewise positive correla- 
tions (P < 0.01) for Cu and K or Mg and Mn were observed. 
It was surprising that Fe was not correlated with P or Zn in the 
cotyledon. 

Meanwhile, for seed coat tissues Fe was positively correlated 
with P and S (at P < 0.01 and P < 0.001, respectively) but not 
with Zn. Potassium and Cu levels were positively and highly 
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FIGURE 1 I Population frequency distributions for iron (Fe) and zinc 
(Zn) concentration in seed coats and cotyledonary tissues in the 
(Cerinza x (Cerinza x (Cerinza x G10022))) advanced backcross 
population as determined by ICP analysis. Arrows indicate phenotypic 



value of recurrent parent Cerinza (C) and the wild donor parent G10022 
(W). The mineral concentration in parts per million (ppm) is found on the 
X-axis, while the number (no.) of recombinant inbred lines (RILs) is found 
on the y-axis. 



correlated (P < 0.001) with those of P and S but less so with 
Mn (P < 0.01). Highly significant negative associations esdsted 
between Fe and B as well as between Mn and Ca in the seed coat. 
Ca was also negatively correlated with K and Mg in seed coat but 
at lower significance {P = 0.01 andP = 0.05, respectively). 

INHERITANCE OF SEED COAT AND COTYLEDONARY IRON AND ZINC 

Using only the Fe and Zn datasets, a total of 7 QTL were found 
with CIM analysis and permutations for significance threshold 
(Table 4 and Figure 2) and 49 for significance at P = 0.05 with 
SPA analysis (Table 4 and Figure 2). Among the QTL identified 
with CIM all were for seed coat Fe and Zn and none were found 
for cotyledonary Fe and Zn. The seed coat CIM — QTL were 
divided into six Zn loci on linkage groups BOl, B02, B07, and Bll 
and one Fe locus on linkage group B04 based on linkage group 
assignments. LR values ranged from 9.49 for the Fe QTL to 29.53 
for one of the Zn QTL. It was notable that the Fe QTL and three 
of the Zn QTL, including one with the second highest LR value, 
were derived from the high mineral wild donor parent, G 10022. 
Three other QTL for Zn concentration, including one with the 



highest LR values and one with the third highest LR value were 
from the cultivated recurrent parent, Cerinza. Two linkage groups 
contained more than one QTL for Zn concentration but in each 
case the two QTL were derived from the same parental allele. For 
example on linkage group BOl, both QTL had positive alleles from 
G10022; while on B02, both QTL had positive alleles from Cerinza. 
The additive effect of the Fe QTL showed that this represents a 
major gene that can provide up to 110 ppm increased seed coat 
mineral concentration, which would change the balance and total 
concentration of Fe in the whole seed. The additive effects of the Zn 
QTL were smaller, ranging from 3.6 to 12.1 ppm increases which 
are also important in total seed mineral concentration considering 
the range of variability for Fe (243 ppm) and Zn (32.43) in seed 
coat. For example the additive effect of 12.1 ppm for Zn is 44.53% 
increase and very similar to that observed for Fe (45.27%). 

Although many SPA — QTL were identified for iron and zinc 
(Table 5) these results agreed with the results from CIM — QTL 
analysis in terms of major locations and source of high iron 
or zinc alleles, validating the previous analysis. SPA results for 
minerals other than Fe and Zn are shown in Table S2. Of the 
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biofortification SPA - QTL, a total of 16 were for Fe in the seed 
coat and 10 were for Fe in the cotyledon. All of the seed coat Fe 
QTL were the result of alleles from G10022, the very high seed 
coat iron parent, and four of the associated markers with highest 
additive effect were on B04. The phenotypic variance explained 
by each of the markers ranged from 7.2 to 25.3 % with the 
marker Pv-at03 and the surrounding markers on linkage group 
B04 being the most important, confirming the results of a major 
locus detected with CIM analysis. The additive effects of substitu- 
tion of the G 10022 allele for any of these markers was an amazing 
1 13 to 171 ppm of iron with BMd9 the most significant marker in 
terms of additivity. 

Other slightly less important SPA-QTL seed coat Fe were 
detected on B02 (at ATA16) and B03 (at ATA26), but these still 
had major additive effects of up to 172 and 100 ppm iron con- 
centration, respectively, for the best marker in each region. The 
phenotypic effects explained by these regions averaged 10 and 9%, 
respectively. A more minor QTL was found on BOl with smaller 
additive effects but still contributing over 55 ppm which is more 
than the average seed concentration of most common beans. The 
phenotypic variances explained by these markers (PV54, 107, and 
233) were around 10%. 

Iron in the cotyledon presented a total slightly smaller num- 
ber of QTL (10) than seed coat QTL (16) with the SPA method. 
However, in contrast to seed coat QTL, most of these SPA-QTL for 
cotyledonary Fe (8) were derived from Cerinza and only 2 were 
from G10022. The position of the QTL for cotyledon Fe were dif- 
ferent than those for the seed coat Fe on some linkage groups for 
example B07 and BIO. On the other hand seed coat and cotyle- 
don QTL overlapped on some linkage groups (BOl, B02, BOS, 
B04, and B08). On linkage group B04, in both cases QTL were 
linked to BMdl6 marker for seed coat and cotyledon, although 
from a different source (Cerinza instead of G 10022, respectively). 
The effects of these QTL were significant as seen in the additivity 
values of 6.6 to 21.9 ppm (up to 56.9% more iron) considering 
the range of variability for cotyledon Fe (38.5 ppm). The QTL 
explained 6.9-15.1% of phenotypic variance. 

In terms of Zn QTL, the SPA method detected 10 QTL for 
seed coat concentration and 12 QTL for cotyledon concentration. 
Among the first of these QTL the associated markers were dis- 
tributed across linkage groups B02, B03, B05, B06, B08, BIO, and 
Bll with most of the high Zn alleles (9) coming from G 10022, 
except for the first one on B02. Their importance however in 
terms of additivity was more minor (6.4-12.7 ppm) although the 
explained variance ranged from 7.3 to a high of 24.6%. In the case 
of the cotyledonary zinc SPA-QTL, 4 were derived from Cerinza 
and 8 from G 10022 and their additive effect ranged from 2.3 to 
9.0 ppm zinc and were found on BOl, B03, B04, B06, B09, and 
Bll. In terms of overlap, the Zn SPA-QTL on linkage group B02 
were different than the CIM-QTL identified on this same linkage 
group. In summary some seed coat Fe or Zn were in similar loca- 
tions as cotyledonary QTL while others were in different locations 
than for cotyledonary QTL. 
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DISCUSSION 

The major results of this study present a quandary for 
biofortification breeding: namely, that both the distribution 
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Table 3 | Pearsons correlation coefficients and significance for ten minerals^ evaluated by inductively coupled plasma (ICP) analysis for seed 
coat and cotyledonary tissue in the advanced backcross population of common beans. 



Coat 


B 


Ca 


Cu 


Fe 


K 


Mg 


Mn P S 


ua 

Pi I 


— U. lUolD 


u.uuyo 












Fe 


P-0.4808'*1| 


0.1085"" 


0.0968"" 










K 


0.2303* 


-0.2054* 


0.2173* 


0.0011"" 








Mg 


-0.1274"" 


-0.2740** 


0.0909"" 


-0.0663"" 


-0.2546** 






Mn 


-0.0133"" 


-0.2959*** 


0.2196* 


0.1159"" 


0.2208* 


0.1239"" 




P 


-0.2396** 


0.0315"" 


0.3075*** 


0.7342*** 


0.4418*** 


0.0613"" 


0.1758"" 


S 


0.0250"" 


0.0613"" 


0.5226*** 


0.2927** 


0.3577*** 


0.1733"" 


0.1404"" 0.6848***^ 


Zn 


-0.0487"" 


-0.0799"" 


0.0559"" 


0.1538"" 


-0.2975** 


0.1154"" 


0.4137*** -0.0767"" 0.0152"" 



Coty B Ca Cu Fe K Mg Mn 



Ca 


0.1024"" 














Cu 


-0.1461"" 


-0.1233"" 












Fe 


0.1192"" 


^ 0.3660 


0.1424"" 










K 


0.0615"" 


-0.1990* 


0.2513*" 


-0.0915"" 








Mg 


0.1419"" 


-0.0495"" 


0.0254"" 


-0.2811** 


0.1431"" 






Mn 


0.0971"" 


0.0631"" 


0.2243* 


-0.0924"" 


0.3531*** 


0.2431*" 




P 


-0.2117* 


-0.1024"" 


0.6508*** 


0.1367"" 


0.3421*** 


-0.0848"" 


0.4886*** 


S 


-0.2208* 


0.0595"" 


0.5806*** 


0.1610"" 


0.0948"" 


-0.0020"" 


0.5513*** 


Zn 


-0.1665"" 


0.1791"" 


0.5358*** 


0.1746"" 


0.1425"" 


-0.0095"" 


0.4922*** 



0.7235*** 0.8094*** 



J 



^Mineral abbreviations for correlation coefficients: B, boron; Ca, calcium; Cu, copper; Fe, iron; K, potassium; Mg, magnesium; Mn, manganese; /? phosphorus; S, 
sulfur; Zn, zinc. 

Tissue Abbreviations: coat, seed coat; coty, cotyledon (and embryo). Significance codes: = 0.001, ** = 0.01, ' = 0.05, ns, not significant, as indicated by dark 
to light shaded boxes, with no significance non-shaded. 



Table 4 | Quantitative trait loci for seed coat iron and zinc concentration identified by composite interval (CIM) mapping analysis in the 
advanced backcross population derived from the wild donor parent (G10022) and the recurrent cultivated parent (Cerinza). 



Trait QTL name 


Chr 


Position 


LR1 


Additivity 


Source 


R2 


Nearest marker^ 






Fe 


4 


0.0901 


9.49 


109.8173 


G10022 


0.087679 


PV-gaatI 


Zn 


1 


0.1701 


27.93 


5.7301 


G10022 


0.062669 


PV139 


Zn 


1 


1.3701 


12.41 


3.6300 


G10022 


0.035964 


PV54 


Zn 


2 


2.4201 


22.69 


4.9477 


Cerinza 


0.042576 


ATA133 


Zn 


2 


3.3001 


29.53 


5.4577 


Cerinza 


0.042892 


PV78 


Zn 


7 


2.4801 


13.95 


12.1286 


Cerinza 


0.274638 


PV35 


Zn 


11 


0.6001 


13.83 


5.4419 


G10022 


0.077907 


BMd33 






Fe None found 


0 


NA 


NA 


NA 


NA 


NA 


NA 


Zn None found 


0 


NA 


NA 


NA 


NA 


NA 


NA 



' LR, likelihood ratio test statistic for Ho:H-i where Ho is the hypothesis of no QTL effect at test position and Hj is the hypothesis of a QTL effect at the test position; 
R2, proportion of variance explained by the QTL at test site; Values have significance at 0.5% probability after 1000-fold permutation tests. 
^The nearest marker is the marker closest to the peak LR score. 



and the inheritance of micronutrients is different in the 
maternally derived seed coat tissue versus in the cotyledonary tis- 
sues which are part of the embryo which will germinate in the 
next sporophytic generation. These results confirm the variability 



in seed coat Fe found by Ariza-Nieto et al. (2007) and Moraghan 
et al. (2002): however, our genetic results are new and differ from 
previous inheritance studies that evaluated micronutrient con- 
centrations in whole seed in other populations (Blair et al., 2009, 
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2010a, 2011) or with the same population (Blair and Izquierdo, 
2012). When dealing with maternally-derived seed coat tissue 
we must think about the seed-producing generation in terms 
of genetics, plant growth and seed development. For any seed 
trait the previous season's growing conditions are important 
too. 



The importance of differences in micronutrient distribution 
in seed coat and the rest of the seed resides in the interaction 
of minerals between each other and with secondary metabo- 
lites such as phytates or tannins that vary in concentration as 
well as between seed coat and cotyledonary tissues (Moraghan, 
2004; Ariza-Nieto et al., 2007). In terms of inheritance of the 
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FIGURE 2 I Continued 
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FIGURE 2 I Location of quantitative trait loci (QTL) for iron (Fe) and zinc 
(Zn) concentration in seed coat. Those QTL identified by CIM analysis in 
the (Cerinza x (Cerinza x (Cerinza x G10033))) advanced backcross 
population are indicated by thick bars acompanied by a QTL name as 



indicated in Table 4. QTL identified for Fe and Zn with SPA analysis as 
indicated in Table 5 are shown as left pointing arrow heads pointing towards 
the most significant markers. Abbreviations for QTL are based on Fe, iron; Zn, 
zinc in seed coat; coat or cotyledon, coty. 
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Table 5 | Quantitative trait loci for seed coat or cotylendon iron (Fe) and zinc (Zn) concentration identified by single point analysis (SPA) in the 
advanced backcross population derived from the wild donor parent (G10022) and the recurrent cultivated parent (Cerinza). 



LG 


Marker 




Significance 




Additivity 


Source 






F 




P 










1 


BM218 


5.23 


0.083 


0.026 


79.16 


G10022 


1 


PV107 


5.63 


0.094 


0.021 


57.80 


G10022 


1 


PV54 


5.06 


0.090 


0.029 


59.94 


G10022 


1 


PV233 


6.81 


0.118 


0.012 


81.59 


G10022 


2 


BMc176 


7.51 


0.128 


0.008 


82.73 


G 10022 


2 


ATA16 


5.55 


0.090 


0.022 


172.43 


G 10022 


3 


PV87 


6.19 


0.100 


0.016 


54.30 


G10022 


3 


ATA 26 


5.05 


0.101 


0.029 


100.33 


G10022 


3 


BM159 


4.12 


0.070 


0.047 


58.88 


G10022 


3 


BM98 


4.96 


0.083 


0.030 


54.60 


G10022 


4 


Pv-at3 


17.36 


0.258 


0.000 


127.50 


G 10022 


4 


PV182 


10.59 


0.172 


0.002 


113.25 


G10022 


4 


BMdie 


6.32 


0.098 


0.015 


132.87 


G10022 


4 


BMd9 


5.26 


0.083 


0.025 


171.37 


G10022 


8 


PV53 


4.17 


0.072 


0.046 


73.31 


G10022 


9 


ATA ME1 


5.03 


0.093 


0.029 


74.97 


G10022 
















Fe COTYLEDON 














1 


PVatccI 


4.74 


0.076 


0.033 


6.61 


Cerinza 


1 


BMc313 


6.15 


0.116 


0.017 


8.54 


Cerinza 


2 


ATA7 


4.22 


0.069 


0.044 


15.78 


Cerinza 


2 


GAT91 


6.36 


0.104 


0.015 


9.58 


G10022 


3 


BMclSO 


4.02 


0.072 


0.050 


7.15 


G10022 


4 


PV-ag4 


4.95 


0.084 


0.030 


7.85 


Cerinza 


4 


BMd16 


5.27 


0.083 


0.025 


12.42 


Cerinza 


7 


BM183 


8.82 


0.143 


0.004 


14.44 


Cerinza 


8 


ATA 289 


7.64 


0.128 


0.008 


21.08 


Cerinza 


10 


BMc159 


9.23 


0.151 


0.004 


21.89 


Cerinza 






2 


BM142 


5.03 


0.088 


0.029 


7.43 


Cerinza 


3 


PV131 


4.57 


0.089 


0.038 


6.40 


G 10022 


3 


BMc333 


4.98 


0.092 


0.030 


10.84 


G 10022 


3 


BMc320 


4.74 


0.079 


0.034 


10.34 


G 10022 


5 


BMd20 


4.59 


0.073 


0.036 


10.34 


G10022 


5 


PV93 


5.60 


0.096 


0.021 


8.44 


G10022 


6 


ATA10 


7.95 


0.124 


0.007 


10.62 


G10022 


8 


ATA 247 


13.07 


0.246 


0.001 


12.73 


G 10022 


10 


BMc9 


4.26 


0.074 


0.044 


6.75 


G 10022 


11 


Bl\/lc368 


4.98 


0.082 


0.030 


17.09 


G10022 
















Zn COTYLEDON 














1 


PV133 


5.53 


0.091 


0.022 


3.12 


Cerinza 


3 


ATA 26 


15.75 


0.259 


0.000 


5.66 


G10022 


3 


BM159 


4.93 


0.082 


0.030 


2.48 


G 10022 


3 


BM98 


6.33 


0.103 


0.015 


2.29 


G 10022 


4 


BMd9 


11.36 


0.164 


0.001 


9.03 


G 10022 


4 


PV182 


9.94 


0.163 


0.003 


4.43 


G 10022 


4 


ATA 143 


8.10 


0.135 


0.006 


3.47 


G10022 


4 


BM171 


5.85 


0.098 


0.019 


4.93 


Cerinza 


6 


BMc371 


4.94 


0.088 


0.031 


2.96 


Cerinza 


9 


ATA ME 1 


10.37 


0.175 


0.002 


3.58 


G 10022 


9 


BM114 


4.22 


0.081 


0.045 


2.51 


G 10022 


11 


BMd33 


4.64 


0.077 


0.035 


3.17 


Cerinza 
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micronutrient accumulation traits, the importance of a difference 
in seed coat and cotyledon resides in the presumably differ- 
ent genes involved in each tissue's mineral accumulation. It was 
not surprising therefore that we found a different set of QTL 
for seed coat Fe and Zn in this study compared to our previ- 
ous study working with whole seed (Blair and Izquierdo, 2012), 
especially in terms of the major CIM-QTL. Seed coat mineral 
accumulation is important to consider for breeding of micronu- 
trient rich beans due to the bioavailability of micronutrients of 
this portion of the seed (Ariza-Nieto et al., 2007) and the GxE 
effects that are found to affect mineral accumulation (Blair et al, 
2010b). 

The genetic results described above showed that most min- 
erals were normally distributed both in seed coat and cotyle- 
don or at least in cotyledon across the population. These data 
proved the multi-genic or quantitative inheritance of min- 
eral concentrations. One exception was the set of values for 
Fe in seed coat, which had a bimodal distribution indicat- 
ing a major CIM-QTL and qualitative inheritance which was 
found to be on linkage group B04. This was in contrast to 
seed coat Zn concentration which was normally distributed and 
had six CIM-QTL, which were found on four different linkage 
groups. SPA-QTL were more abundant and will be discussed 
as well. 

Blair and Izquierdo (2012), with the same population but 
analysis of whole seed, found the same CIM based QTL 
for Zn accumulation on B07 near marker PV35 but dif- 
ferent QTL for Fe concentration or content on B07 and 
BOS, respectively. Some SPA based QTL for Fe and for Zn 
from Blair and Izquierdo (2012) overlapped with the seed 
coat or cotyledon based SPA - QTL on BOl, B04, BIO, 
andBll. 

The QTL for zinc concentration on Bll were especially impor- 
tant in another Andean x Mesoamerican cross study (Blair et al, 
2009). The QTL for zinc concentration on BOl near PV139 
may have been influenced by the fin locus as was found by 
Cichy et al. (2009) for another cross involving a determinate 
Andean bean, similar in growth habit to Cerinza. The fin locus 
was shown to control the determinacy but not the height of 
bean plant growing shoots (Chavarro and Blair, 2010) and in 
the advanced backcross population analyzed here was observed 
to affect total biomass production and plant size (Izquierdo 
et al., unpublished results). One possible hypothesis could be 
that larger indeterminate Fin Fin plants accumulated a large 
amount of zinc in vegetative tissue and this zinc was avail- 
able for translocation to the seed. Meanwhile, the shorter and 
smaller biomass fin fin plants that are determinate in growth 
habit would accumulate less zinc in both vegetative and repro- 
ductive tissues. A different Zn QTL found on linkage group 
B02 near PV78 and derived from Cerinza was in close prox- 
imity to the QTL ZnPoAAS2.1 found by Blair et al. (2010a) in 
a Mesoamerican x Mesoamerican cross and another QTL Zn- 
AAS2c found by Blair et al. (2011) in an Andean x Andean 
cross. 

On the other hand, many additional SPA-QTL from this study 
appear to be novel based on their evaluation in the different tis- 
sues. For example, the seed coat Fe QTL on B04 derived from 



G 10022 was not detected near any Fe concentration QTL in the 
cultivated Andean x Andean or Mesoamerican x Mesoamerican 
populations studied by (Blair et al., 2010a, 2011) and may 
be specific to wild bean sources. This major QTL from link- 
age group B04 probably influenced the observation of binomial 
population distribution in the population for seed coat iron 
concentration. 

The lack of cotyledonary CIM-QTL for iron or zinc in the 
present study may have been a reflection of the smaller differences 
between parents for the Fe and Zn concentrations in this tissue, 
although SPA-QTL analysis identified a good number of markers 
with significant effects. This may show that epistatic interactions 
are sometimes important in identification of QTL for the min- 
erals due to the dependence of iron translocation to the seed 
on the amount of iron uptake by the roots into the plant. More 
QTLs were found with the SPA method due to the lower proba- 
bility threshold (P < 0.05) used in that analysis compared to the 
CIM method where thresholds were determined with permuta- 
tions. The lack of cotyledonary CIM iron QTL was a result of 
high thresholds found for this variable's analysis. Despite this, the 
cotyledonary QTL where in similar locations to QTL found in 
the whole seed for the same population by Blair and Izquierdo 
(2012). 

From a breeding perspective, the lines with high seed coat 
Fe concentration might be of interest if this Fe is shown 
to be bioavailable although initial results suggest that seed 
coat iron especially in colored beans is not very bioavail- 
able (Ariza-Nieto et al., 2007). It was notable that some 
lines of the population had seed coat Fe concentrations of 
more than 250 ppm, while seed coat Zn concentrations above 
50ppm might be of interest. The microsatellite that were 
associated with QTL for seed coat mineral concentration 
could be used for marker assisted selection (MAS) of these 
traits as was suggested for whole seed concentration QTL by 
Blair and Izquierdo (2012). 

Considering that the seed coat makes up approximately 10% 
of the seed the amount of iron in the seed coat can have a 
large impact on the amount of iron in the entire seed. In pop- 
ulations like the advanced backcross lines, the amount of seed 
coat is fairly uniform but the amount of seed coat Fe was not. 
In a typical recombinant inbred line population that segregates 
for seed size, the percentage seed coat, the ratio of seed coat to 
total seed and the amount of seed coat Fe would all be vari- 
ables to study. The result of the high iron concentration in the 
seed coat and the high percentage of seed coat to seed weight 
is that seed coat Fe at high concentrations affects the over- 
all average Fe concentration of the seed quite substantially. For 
example, from this study and our previous results (Blair and 
Izquierdo, 2012), the total amount of seed iron can increase 
or decrease by up to 40 ppm based mosfly on seed coat Fe 
content. 

The results we observed for correlation among minerals for 
cotyledons are different than Fe-Zn relationships observed pre- 
viously for whole seed in several populations (Blair et al., 2009, 
2010a, 2011) and in the same advanced backcross population 
(Blair and Izquierdo, 2012). Meanwhile seed coat results showed 
interactions of Fe and P which are new information for seed coat. 
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but there was no correlation for Fe and Zn, in contrast to Blair 
et al. (2009). A relationship of Fe and P in seed coat was surprising 
considering that phytate which binds Fe should not be found in 
high concentrations in the seed coat but rather in the cotyledons 
(Ariza-Nieto et al, 2007). 

The Fe and P correlation in seed coat but not in cotyledons 
can be influenced by the fact that concentrations of P in the 
seed coat are only one-tenth that of P concentrations in the 
cotyledons and that in seed coat P may be binding with tan- 
nins and other seed coat substances that also bind Fe (Blair 
et al., 2012). The lack of correlation in the cotyledon between 
P and Fe might be because Fe in legumes is targeted to vas- 
cular cells where starch also accumulates in amyloplast and is 
found bound to ferritin (Cvitanich et al., 2010). Perhaps phy- 
tates are not important in sequestering iron in the seed but rather 
only as anti-nutrients in the human digestive system (Ariza- 
Nieto et al., 2007). On the other hand P fertilization is known 
to affect Zn uptake by plants (Moraghan and Grafton, 2002) 
and Zn QTL detection (Cichy et al., 2009). The number of Zn 
QTL contrasts with single gene inheritance found by Cichy et al. 
(2005) in navy beans, but agrees with multiple gene inheritance 
found in later studies (Blair et al., 2009, 2010a, 2011; Cichy et al, 
2009). 

Our hypothesis of a major gene for Fe in seed coat on B04 
is the first time a single gene has been proposed for control of 
this mineral's accumulation in the seed coat of common bean and 
therefore might be a target for gene cloning and characterization. 
It would be interesting to know if the seed coat accumulation 
of Fe is based on a seed coat expressed gene or a gene that 
diminishes loading of iron into the seed's embryo/cotyledons. 
Based on co-localization or the lack thereof, ferritin is unlikely 
to be the mechanism of iron accumulation in seed coat but 
binding with seed coat tannins could be possible (Ariza-Nieto 
et al, 2007; Cvitanich et al, 2010). The location of the QTL on 
the long arm of B04 is interesting as nearby genes for phyto- 
hemaglutinnins are expressed only in seeds (Blair et al., 2010c). 
The clustering of genes expressed in the same tissue is often 
typical. 

Whatever the case, the specificity of the gene to the wild 
bean source is interesting and could be relevant for gene dis- 
covery in other legumes, such as the model species Medicago 
truncatula, where transcriptome analysis has shown the large 
number of genes (over 30,000) expressed in the seed coat (Verdier 
et al., 2013). Before this study, the high iron of wild beans could 
be questioned as an artifact based on the high ratio of seed 
coat to total seed size in these very small seeded seeds (less 
than 10 g per 100 seed). If this is the case it would be diffi- 
cult to use wild beans as a source of high iron. However, it 
appears that wild beans might preferentially accumulate iron 
in their seed coats compared to cultivars. Alternatively, cul- 
tivars might have been selected to accumulate less iron in 
their seed coats than wild beans as part of domestication and 
development of large seed sized domesticates. Selection pro- 
cesses may have been different in large-seeded Andean vs. 
small-seeded Mesoamerican beans or in very thick-seed coated 
species such as scarlet runner bean and year-long bean (Singh, 
2001). 



The idea that human selection for these invisible micronu- 
trient traits has been active is intriguing. We could postulate 
that this selection was conducted through proxy mechanisms 
such as seed coat color, palatability, cooking time or notice- 
able health benefits. Seed coat color is an obvious trait that 
differs greatly in common bean cultivars with some of these col- 
ors associated with tannin or anthocyanin accumulation (Caldas 
and Blair, 2009; Diaz et al, 2010). Perhaps humans instinc- 
tively were attracted to low tannin varieties of white and yel- 
low beans in some cultures or developed methods to remove 
tannin in red and black beans by soaking, thus also affecting 
palatability. 

Palatability is hard to measure, but cooking time parame- 
ters have been evaluated in various studies that also have looked 
at micronutrient concentration. For example, Saha et al. (2009) 
found that zinc concentration was negatively correlated with 
swelling and hydration capacity as components of cooking qual- 
ity. A few high iron genotypes from that study were also hard 
to cook. The amount of Fe in the seed coat might have impli- 
cations for the cooking method recommended when aiming 
to preserve the nutritional quality of common beans (Ariza- 
Nieto et al., 2007; Carvalho et al., 2012). As discussed, some of 
the QTL specific to the seed coat concentration for Fe and Zn 
were the same as those identified for overall mineral content 
and might lead Fe to be more likely to stay in the seed when 
cooked. 

The idea that farmer-consumers selected for micronutrient 
density in the parts of the seed that were more bioavailable in 
their cultivars compared to their wild relatives is not far-fetched: 
iron consumption may be noticeable in a feeling of more energy 
from avoiding anemic status which would have been connected 
to the foods eaten. The diet of Amerindians from the tropi- 
cal New World who domesticated common beans was unlike 
the diets of Indigenous people to the north or south which 
were very low in animal protein sources. Therefore, the tra- 
ditional diet of pre-Colombian Middle American and Andean 
cultures are almost exclusive in vegetable protein sources with 
basic staples such as maize, beans and squash that were bet- 
ter in amino acids, micronutrients and vitamins than the sta- 
ples of "Old World" consumers, who had carbohydrate-rich 
millets, sorghum, rice and wheat. Beans are known as the 
"meat of the poor." Perhaps they should be better known as 
"meat of the Americas" and as "a substitute for animal pro- 
tein" as many poor people or vegetarians in Africa, Asia and the 
Americas are well aware of, based on their high consumption of 
beans. 

In conclusion, further molecular and physiological approaches 
are needed to explain the mechanism for iron accumulation in 
the seed coat versus the cotyledons but this study provides a 
starting point for gene and tissue analysis of seed coat minerals. 
Additional work should look into why there is a very different 
genetic control of seed iron and zinc tissue distribution in wild 
vs. cultivated beans, with higher concentrations of iron accumu- 
lating in the seed coat of wUd beans, but lower concentrations 
in their cotyledons, and the reverse being the case for culti- 
vated beans. In addition to a major gene on B04 for seed coat 
Fe in both CIM and SPA analysis, it was notable that the most 
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important of the seed coat Zn SPA-QTL was from a region 
of the genome on linkage group Bll, that has been impor- 
tant in other studies (Blair et al, 2009; Cichy et al., 2009) 
and that contains ZIP type transporters, which have a role in 
the uptake and transport of iron in plants (K. Cichy, pers. 
Communic). Some other regions on BOl, B02, B03, B04, B09, 
and BIO overlap for QTL with those of studies of Fe and Zn 
accumulation in whole beans (Blair et al., 2009, 2010a, 2011; 
Cichy et al., 2009). The distribution of QTL for other miner- 
als or the effect of seed weight QTL on B02, B03, and B09 
might be worth looking at in more detail in respect to Fe 
and Zn accumulation. The implications of all this work for 
biofortification of beans should all be considered in molecu- 
lar breeding and physiological analysis. Genetic studies along 
with candidate gene analysis provide the framework to under- 
stand the physiology of iron uptake into the plants, trans- 
fer to the pod and seed and accumulation in different seed 
tissues. 
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